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Abstract 

A nonsingular localized static classical solution is constructed for standard Einstein gravity coupled 
to an SO(3) x SO(3) chiral model of scalars [Skyrme model]. The construction proceeds in three 
steps. First, an Ansatz is presented for a solution with nontrivial topology of the spacetime 
manifold. Second, an exact vacuum solution of the reduced field equations is obtained. Third, 
matter fields are included and a particular numerical solution is found. This numerical solution 
has a negative effective mass, meaning that the gravitational force on a distant point mass is 
repulsive. Most likely, the origin of the negative effective mass traces back to the "surgery" needed 
to create the "defect" from Minkowski spacetime, but this process involves topology change and 
lies outside the realm of classical Einstein gravity. 
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I. INTRODUCTION 

It can be argued on general grounds that the small-length-scale structure of quan- 
tum spacetime is nontrivial. This structure has been called a quantum spacetime foam. 
Over larger length-scales (lower energy-scales), an effective classical spacetime manifold 
emerges and the crucial question is if that effective spacetime is perfectly smooth or not. 
Particle/ wave propagation over Swiss-cheese-type manifolds has been studied and the prob- 
lem can, in principle, be solved exactly 

aa 

. The simplest such Swiss-cheese-type manifold 
has identical static "defects" (alternatively called "holes" or "knots"), where each defect 
provides nontrivial topology. The particular defect considered has, however, a divergent 
(delta-function-type) Ricci curvature scalar and does not solve the vacuum Einstein equa- 
tion. 

The goal, now, is to construct a nonsingular defect solution by use of appropriate matter 
fields. For the defect topology at hand (holes with antipodal points identified; see below), it 
has been suggested ^] to use the gravitating SO (3) Skyrme model [g, 7] with an appropriate 
interaction term [8] allowing for negative energy-density contributions. The present article 
uses a new Ansatz for the fields and gives special attention to the behavior of the reduced 
field equations at the defect core. We are, then, able to obtain a nonsingular defect solution 
in the gravitating 50(3) Skyrme model. 

Remarkably, the obtained solution displays the phenomenon of antigravity, i.e., a nega- 
tive effective mass. As will become clear later on, precisely this antigravity property and 
the presence of negative energy densities allow a particular singularity theorem to be 
circumvented. 

II. MANIFOLD 

The 4-dimensional spacetime manifold considered in this article is static for appropriate 
coordinates, 

I 4 = lxM 3 . (2.1) 

The nontrivial topology appears in the 3-space M 3 , which is, in fact, a noncompact, ori- 
entable, nonsimply-connected manifold without boundary. Up to a point, M 3 is homeomor- 
phic to the 3-dimensional real-projective plane, 

M 3 ~ RP 3 - {point} . (2.2) 

Further details can be found in Refs. 

aa 

. Here, only the necessary information is given. 
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Figure 1: Three-space M 3 obtained by surgery from M 3 : interior of ball with radius b removed and 
antipodal points on boundary of ball identified (as indicated by open and filled circles). 



For the explicit construction of M 3 , we perform some local surgery on the 3-dimensional 
Euclidean space E 3 = K 3 . We use the standard Cartesian and spherical coordinates on M 3 , 

x = |x| x — (x , x 2 , x 3 ) = (r sm9 cos0, r sin # sin 0, rcos^) , (2-3) 

with x rn G (-00, +00), r > 0, 9 G [0, 7r], and G [0, 2tt). Now, M 3 is obtained from M 3 by 
the removal of the interior of the ball Bj, with radius b and the identification of antipodal 
points on the boundary S b = dB b . With point reflection denoted by P(x) = —x, the 3-space 
M 3 is given by 

M 3 = {x G R 3 I (|x I > 6 > 0) A (P(f) S f for |f| = &)} , (2.4) 

where = stands for point-wise identification (Fig. [Q. 

The single set of coordinates (I2.3P does not suffice for an appropriate description of M 3 . 
The reason is simply that two different values of these coordinates may correspond to a 
single point. For example, x = (b, 0, 0) and x = (—6, 0, 0) correspond to really the same 
point of M 3 . A relatively simple covering of M 3 uses three sets of coordinates (also called 
charts or patches), labeled by n = 1, 2, 3. Each coordinate patch surrounds one of the three 
Cartesian coordinate axes. The sets of coordinates are denoted 

(X n ,Y n ,Z n ), for n = l, 2, 3, (2.5) 

which are, despite appearances, really non-Cartesian coordinates. Specifically, the set of 
coordinates surrounding the x 2 -axis segment with \x 2 \ >b is given by 

X 2 = < > fOT ° <0<7r ' (2.6a) 

!> — 7r for 7r < (f) < 2ir , 

. . • — b for < 6 < tc , , . 

Y 2 = { * (2.6b) 

b — r for 7r < <p < 2-7T , 

e o<*<,, (2&) 

7T — 6 for 7T < < 27T , 
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with ranges 

X 2 e (0, tt) , (2.7a) 

Y 2 e {-oo, oo) , (2.7b) 

Z 2 e (0, tt) . (2.7c) 

The other two sets, (Xi, Yx, Z\) and (X 3 , Y" 3 , Z 3 ), are defined similarly. 

In the following, we consider spherically symmetric fields and it suffices to consider one 
set of coordinates, which we take to be n = 2 set (12. 6p . The notation is simplified as follows: 

(X, Y, Z, T) = (X 2 , Y 2 , Z 2 , t) , (2.8) 

where the time coordinate has been added in order to describe the spacetime manifold M4. 

III. FIELDS AND ACTION 

The spacetime manifold (12. ip of the previous section is now equipped with a metric g^, 
whose dynamics is taken to be governed by the standard Einstein-Hilbert action jloj |. In 
addition, there is assumed to be a scalar field fl e SO (3), with self-interactions governed by 
a Skyrme term in the action p] and by another term [8| whose coupling constant 7 is taken 
to be positive, allowing for negative energy-density contributions. 

Specifically, the combined action of the pure-gravity sector, labeled 'grav,' and the matter 
sector, labeled 'mat,' is given by (c = h = 1) 

S = I d X y/ g I £g r av, EH 4~ -£-mat,kin 4~ -^rnat, Skyrme 4" £mat, metastab ) ; (3.1a) 

where 

£ ^ EH = 16^*' (3 ' lb) 
f 2 

^ mat ,kin = — trfu^a/") , (3.1c) 

£ mat , Skyrme = tr( [u^ U„] U v ] ) , (3. Id) 

1 / \ 2 

Anat, metastab = \ tl (u^U^) J , (3.1e) 

in terms of an so(3)-Lie-algebra valued one-form u = Co> M dX^ given by 

w^ = n -1 9 M fi. (3. If) 
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The SO (3) x 5*0(3) global symmetry of the matter sector is realized on the scalar field as 
follows: 

n ->■ s L ■ n ■ Sn 1 , (3.2) 

where Sl, Sr G 50(3) and the central dot denotes matrix multiplication. As mentioned at 
the end of Sec. [TTJ the measure d 4 X in the integral (13.1 al) corresponds to only one of the 
three coordinate charts needed to cover M4. 



IV. ANSATZ AND FIELD EQUATIONS 



Ansatz 



A spherically symmetric Ansatz for the metric is given by the following line element: 

dg 2 = _ e 2v(W) dT 2 + e 2~X(W) dY 2 + ^2 + ^2 z ^ ? ^ ^ 

W = b 2 + Y 2 , (4.1b) 

where, as explained at the end of Sec. [Til we only show the coordinates of one chart with 
Y G (— oo, oo). The scalar field is given by the hedgehog-type Ansatz 

HQ, 



Q = cos [F (r 2 ) ] 1 - sin [F (r 2 ) ] x- S+ (l — cos [F (r 2 ) ]) 



X g) X . 



Fib 2 ) 
Si 



71 . 



/ o o o\ 

1 
0-10 



So = 



0-1 



1 



J 



( o i o\ 


-10 



(4.2a) 
(4.2b) 

(4.2c) 



where, in components, (x <S> x) ab = x a x b . Remark that, because of the boundary condition 
(14.2bj) at \x\ = b, it is possible to use the single coordinate chart ( 12. 3ft with the further 
identification r 2 = b 2 + Y 2 = W for the coordinates used in the metric ( 14. ip . In other 
words, the topology of M3 (see Fig. [TJ is trivially consistent with the hedgehog field having 
boundary value F(b 2 ) = tt. 

The arguments of the tilde-functions in the above Ansatze have dimension length-square. 
Related functions with lengths as arguments can be defined as follows: 



(V& 2 + F 2 ) = v (b 2 + Y 2 ) 



(4.3) 



and similarly for A and F . The functions without tilde resemble those of the previous 
literature [s-8], but we prefer to work with the tilde-functions. 
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The relevant nonvanishing components of the Riemann tensor for the Ansatz (14. ip are: 

R T yTY = -2 [u 1 + 2Y 2 {-V k'/k + V + v' 2 )\ , (4.4a) 

R T ZTZ = -2Y 2 u'/k 2 , (4.4b) 

R T xtx = sin 2 ZR T ZTZ , (4.4c) 

R ZYZ ~ (b 2 + Y^ ' (4 - 4d) 

R Y XYX = sin 2 ZR Y ZYZ , (4.4e) 

= sin 2 Z f 1 - m Y *J \ , (4.4f ) 



(b 2 + Y 2 )n\ 

where the prime stands for differentiation with respect to W. The components not shown 
in (14.41) either vanish or can be computed from the ones above by the usual symmetry 
properties. 

B. Reduced field equations 

The derivation of the variational equations is straightforward. 1 Henceforth, we will use 
the following dimensionless model parameters and dimensionless variables: 

rj = 8nr] = 8nG N f 2 , (4.5a) 

w = (ef) 2 W = (y ) 2 + y 2 , (4.5b) 

y = efY, (4.5c) 

Vo = efb. (4.5d) 

Furthermore, the following functions are introduced: 

k(w) = exp [A(w)] , (4.6a) 

Jl(w) = exp . (4.6b) 

The reduced Einstein equations and matter field equation will be given in Appendix |A3 



1 The variational equations for 7 = can also be obtained by making appropriate rescalings (F — s- F/2 and 
e 2 — > e 2 /2) and substitutions (v — > u, A — > A, and F — > F) in the equations of Ref. Q. 
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From these equations, one obtains three ordinary differential equations (ODEs): 

u'(w) = 4^gj - ^ + vQi [F{w), *(«,), ?(«,)] , (4.7a) 

2 ' — ' / \ 2 

k'(w) = k(w) W ± ~ *y + ^Q 2 [FH, *(«;), u(w)] , (4.7b) 
4w (u> — Uq) l j 

F» = Q 3 [F(w),k(w),u(w)}, (4.7c) 

where the prime now stands for differentiation with respect to w and Q n are certain func- 
tional given in Appendix |A] These ODEs are to be solved with the following boundary 
conditions: 

F(0) = F(oo) =tt, (4.8a) 

k(0) = 0, z7(oo) = 0. (4.8b) 

The field configuration (14. 2 p with boundary conditions (14.81) has winding number zero 
[recall ( 12. 2 j) and the further homeomorphism 5*0(3) ~ MP 3 ]. It is, therefore, not a genuine 
Skyrmion [which has F(oo) = and winding number one]. Our matter solution may very 
well turn out to be unstable, but that is not important for one of the main goals of this 
article, i.e., finding a possible origin of antigravity effects (see Sec. IVT) . 

C. Power-series solution 

An analytic solution near the defect core can be obtained as a power series. The first 
terms are as follows: 

F(w) = ix -f(w- yl) + ■■■ , (4.9a) 



k(w) = ky/w -$ + ••• , (4.9b) 

p(w) = u - n (w - yl) H , (4.9c) 

with coefficients /, k, and n, and the arbitrary constant uq. Given the value for k, the 
reduced field equations fix the values (denoted by an overbar) of the other two constants: 



4>yyo 2 



n — k' 



1 , Q ~ 377 9^ 6r? 4r?7 3r]y 2 
-- + 377 + ^ 



(4.10b) 



4 ' 7 ' 1/0 2 iyo 2 yo 2 87 

where 7 is assumed to be positive and not too large, so that 12 — 8 7 + 3 yo 2 > 0. 

With the leading terms of the above power-series solutions, it is easy to verify that the 
components of the Riemann tensor (I4.4p are finite at the defect core y = 0. 



7 



V. EXACT VACUUM SOLUTION 



For vacuum matter fields, 

F(w) = 0, (5.1a) 



the ODEs (I4.7aj) and (I4.7bj) have Qi = Q2 = and can be solved exactly for boundary 
conditions (14.8bj) . 



k(w) = exp [AH] = . (5.1b) 

jj,(w) = exp [^(itf)] = \Jl — 1/ y/w , (5.1c) 

with a dimensionless constant £ and the definition w = y 2 + (yo) 2 in terms of the 'radial' 
coordinate y G (—00, +00) and the defect parameter y > (the corresponding dimensional 
parameter b is also nonzero and positive). The same solution is obtained for the case 
F(w) = 7r and rj = 0. 

At this point, we can already make two general remarks. First, the vanishing of the 
metric component H 2 at Y = implies that the defect of Fig. [T]at a fixed time T has zero 
physical extent along the radial direction. 

Second, the real k and Jl functions from (15.11) . for large enough values of \£\, cover the 
whole of the manifold M 4 only if £ is negative: 

£ 2 >(y ) 2 => £<-y <0- (5.2) 

In turn, this implies having anti- gravity, because asymptotically (w ~ y 2 ) the solution 
behaves as follows: 

k 2 ~ 1/(1-^1), (5.3a) 

Jl 2 ~ (1-^/M), (5.3b) 



which is to be compared to the standard Schwarzschild metric [10[ with positive mass M, 
having the respective components 1/(1 — IG^Mjr) and (1 — IG^Mjr). For the vacuum 
solution ( 15. ip with £ < 0, a point mass far away from the defect core will not be attracted 
towards it but repulsed. 

For £ < yo, all Riemann-curvature-tensor components ( 14. Ah are finite over the whole 
manifold, also at the defect core, y = or w = (yo) 2 . Specifically, we find for the relevant 
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nonvanishing Riemann tensor components: 

R T yty = (e 2 f) 4r7=4v > (5.4a) 

R T ZTZ = (5.4b) 

f 

R z xzx = sin 2 Z^. (5.4d) 



w 

More importantly and for all values of £, the scalar obtained by contraction of the Riemann 
tensor with itself, 

S 2 = R^paR^ = 12 —„ > ( 5 - 4e ) 

W a 

remains finite over the whole of the manifold M 4 , because w > {yo) 2 > 0. This behavior 
contrasts with that of the Schwarzschild metric over R 4 , for which S 2 diverges at the point 
r = 0. 2 

A brief discussion of the geodesies from the metric ( 14. ip with functions ( I5.1bj) and (I5.1cp 



is given in Appendix [Bj A more extensive discussion of geodesies is scheduled for a follow- 
up paper, while another follow-up paper intends to discuss the global structure of the new 
vacuum solution over M 4 . 3 

Note, finally, that the regular flat-spacetime metric is given by (15.1bj) and (15. left with 



£ = 0. The actual value of the free parameter I in the metric from ( 15. ip will have to be 
determined by adding matter fields (the same applies for the determination of M in the 
standard Schwarzschild solution). An interesting question is if it is possible to get negative 
£, corresponding to antigravity. 

VI. NUMERICAL SOLUTION WITH MATTER FIELDS 

We use the following values for the model parameter 7 and the defect parameter y : 
7 = 3/2 + A , A > 0, (6.1a) 



To = V[7(8 + 5 2 )-12]/3. (6.1b) 



2 With < b < £/(ef) = 2G N M, the manifold (J2H and metric ([13]) . (jglbl) . (|5Jcj) provide a "regulariza- 
tion" of the standard Schwarzschild solution. 

3 Electromagnetic wave propagation over this spacetime can be calculated with the methods of Refs. [3, |4|. 
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These parameters for 5 — > are consistent with having zero coefficient / and positive finite 
coefficient k in the series solution ( I4.9p . according to ( I4.10p . The numerical solution of 
the ODEs will be obtained for w G [u> m in, w max ], with w m \ n = (yo) 2 + e 2 and large enough 
w max . The boundary conditions at Wb cs = w min follow from the analytic expressions (I4.9p . 
Anticipating the discussion of Sec. IVIIIt the coupling constant (I4.5a|) is taken infinitesimal, 
1] = e 2 . 

For parameters (I6.ip with 5 = and A = 1/100 and for a particular value of the coefficient 
k, the numerical solution is given in Figs. [2]and|3j The dimensionless energy density is seen 
to be negative at the defect core [with too,def = 4e 2iy ° (l + y 2 — 47/3)/ (|/o) 4 ] and to have a 
positive tail with diverging total energy [see Sec. I Villi for further discussion]. 

From the asymptotic behavior of the functions shown in Fig. [2j the parameter i as defined 
in (15.31) is found to take the value 

£^-28, (6.2) 

also for smaller values of e. The numerical result (I6.2p suggests that it is possible to realize 
the vacuum solution (15. ip with a negative mass parameter £. Still, the numerical analysis 
reported in this article is only exploratory and needs to be confirmed. 



VII. DEFECT MASS 



According to (I5.3p . the parameter £ of the vacuum solution (15.11) is proportional to twice 



the Schwarzschild mass (see Sec. IVIIII for further discussion of the energy sea 



es of the 



model). Another useful definition of the mass is the so-called Komar mass 1CJ, For the 
Ansatz (14. ip . the dimensionless Komar mass at infinity becomes 



^Komar = lim ^Komar(w) = Jim — — 2w \ W - y\ v' (w) . (7.1) 

«i—s>oo w— >oo K[W) v 

This mass comes from an integral over a distant sphere (having radius square w = w), but, 
for the vacuum solution, it can be transformed |5( into an integral over the defect sphere 
(having w = y$) with the result: 

WKomar = "iKomar(l/o) = ■ ( 7 - 2 ) 

Even with matter fields included, as was the case for the numerical solution of Sec. IVIt 
the negative contribution to the Komar-mass from the "defect" at r = b can still dominate 
the standard (positive) contributions from gravity and matter in the bulk. Precisely this net 
antigravity property and the presence of negative energy density at the defect core invalidate 
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Figure 2: Numerical solution of the ODEs (|A4p over the w interval {iu m in> i^max} = {2/75 + 
e 2 , 10 6 } with e = 10~ 4 and model parameters {7, 77, yo} = {151/100, 10~ 8 , ^2/75}. The boundary- 
conditions at Wbcs = Wmin = Vq + <? follow from the analytic expressions (14. 9p with coefficients 
{k, n, /} = {a/8/75, -(1/4 + 37 vr/12500000) 8/75, 0} and constant u adjusted to give ?(w max ) = 
-X(w max ). 
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Figure 3: Dimensionless energy density too = 2oo/(e 2 / 4 ) for the functions from Fig. [2l 

the Singularity Theorem 2.1 of Gannon |[9|, which applies to spacetimes with asymptotically- 
flat nonsimply-connected spacelike slices. This suggests that having a nonsingular defect 
solution may result in a negative mass, also for matter sectors which are different from the 
particular Skyrme model considered here. 

The defect mass discussed up till now is the gravitational mass. It remains to be seen 
what the inertial properties of the defect are. 



11 



VIII. DISCUSSION 



The numerical solution from Sec.|VT]has asymptotically a nonzero energy density, (Fig. [3]) 



Too oc - e 2 / 4 ~ f ~ , (8.1) 
w r z 

where we have used w ~ y 2 ~ (e/r) 2 in terms of the standard radial coordinate r and the 
parameters e and / from the action (13. ip . The asymptotic value of the Komar mass of the 
numerical solution is 

MKomar ~ i/2 , (8.2) 



where the coefficient £ is obtained from the behavior 2A ~ — 2u ~ — £/\y\ in (15. 3p . 
Of interest to physics are, therefore, the following mathematical limits: 

7^0, (8.3a) 
e -»■ oo, (8.3b) 
e / — >■ const > , (8.3c) 

while is kept fixed (recall c = 1). The limit (I8.3aj) makes for vanishing energy density 
(18. ip . whereas the limit (I8.3cp keeps the Komar mass (18. 2 p finite. 

The numerical solution of Fig. [2] has a Schwarzschild-like behavior (I5.3P with i < 0. The 
corresponding nonsingular defect in the limits (18. 3p then produces a finite negative mass, 
i.e., antigravity, while providing for a negligible local energy density (the total scalar- field 
energy up to a fixed physical distance also vanishes). However, as mentioned at the end of 
Sec. Wl\ the numerical analysis of the present article still needs to be verified and extended. 

Note that, strictly speaking, the 5*0(3) scalar field is absent in the Standard Model. But 
the nonlinear sigma model resurfaces if the gauge fields are eliminated and the Higgs-field 
modulus is frozen. Hence, the toy model of matter fields considered in the present article 
is not completely removed from realistic physics. As such, there may even be a connection 
with ideas linking the quantum structure of spacetime (having an energy scale -Epianck 



sjhc 3 /Gm and a length scale frc/Epi^^) to the Higgs-boson and top-quark masses 12]. 
We refrain from making even wilder speculations about the cosmological constant. In fact, 
the main focus of this article remains purely theoretical, namely, finding a nonsingular 
spacetime-defect solution and a possible mechanism for the origin of antigravity effects. 
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Appendix A: ODEs 



In this appendix, the relevant ordinary differential equations (ODEs) are given in detail, 
whereas they were only given in abbreviated form in ( 14. 7p . It turns out to be helpful to 
define the following dimensionless auxiliaries: 



Ml 



1 7 sin — 

3 V I 2 



(Ala) 



u 2 



(Alb) 



«3 = V 



4 \ ~ / 4 \ ~ / 4 4 
1 - - 7 1 cos(2F) - 4(l--7 + wlcosF + 3(l--7+-w 



- 2w . (Ale) 



The reduced Einstein equations can now be written in the form G Y — SttG^T Y , G Y ~ 
G T T = 8ttG n (T y y -T t t ), and G z z = 8ttG n T z z [it is also found that G\ = G z z , T\ = T z z }. 
In terms of dimensionless quantities, these reduced Einstein equations are: 

(w - yl) e~ 2X (Awu' + 1) _ 



4?7 



w 



1 + ^) sin 2 (-)+(«;- Vl) e- 2X F' 2 ( -2u 2 + 1W (w - y 2 ) e~ 2X F 



2(w-yl)(X' + u')-^ 



w 



1 



8^7 (w - yl) F' 2 [ u 2 - pw (w - y 2 ) e~ 2X F' 2 



-e~ 2X [2w{w- y 2 ) ( \> V - v' 2 - v") + (w - y 2 ) X' - (2w - y 2 ) V - 



2w 



(A2a) 



(A2b) 



-VUism- I — 1 -w V (w-yl) e~ 2X F' 2 (l - - 7 (w - yl) e~ 2X F' 2 



(A2c) 
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The reduced matter field equation is 



0=1 



2«i 



1 1 „2A 



W 



e ZA sin F 



-l§\lw-zlF' 



u 2 [(w-zl) - A 1 ) + l) + - (to - zg) 



16 



1 (w- zl) 3 ' 2 e~ 2 ~ x F' 3 \w (w - zl) fa - jA - 3w + z 2 



20 



- 2 (l-^) ( W -, 2 ) (l-4^-z 2 ) F' 2 sinF 



(A3) 



-16 (w - zl) 3 ' 2 F" (u 2 -~f W (w- zl) e- 2X F' : 

Using k(w) as defined in fl4.6al) and solving for u', k', F", one obtains the final ODEs: 

~, k 2 1 riK 2 9 f F\ , , _ ~, 9 

4 (w - yjj) 4w u> (w - j/5) y 2 y 

-7^w (w - yl) k~ 2 F' A , (A4a) 

(A4b) 

o 

1 



K = K 



w + yl — k 2 



24 ( U2 — 7W (w — yl) k 2 F' 2 

m J- v lfn + f ) 8i " * - ^ + 4M2 + 2w " ^ I F ' 



1 - AJw -yl) 1 - -7 sinF F 



-4 7 (m 3 + (-Aw + 3yl) k' 2 ) F' 3 -16jr}w(w- yl) k~ 2 u 2 F' 5 



(A4c) 



Appendix B: Vacuum-solution geodesies 



For the vacuum solution of Sec. |V] the nontrivial Christoffel symbols are 



ef 
1 

~f 
1 



1 rprp 



YY w (v^J - t) yjw - yl ' 



Vw~^li/2 
(y/w~-l) 1/2 



w^Jw- yl 
wy 2 -(w + yl) 1/2 



(Bla) 
(Bib) 

(Blc) 
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The last two Christoffel symbols are seen to diverge at the defect core (w = y^), but 
the particle motion can still be regular, as we will now show by a simple example [using 
dimensionless Christoffel symbols 7^, = (ef)' 1 T^ c }. 

Consider, in fact, the geodesic equation for a particle moving solely in the y-direction and 
denote the nonvanishing dimensionless velocity components u l = dt/dr and u y = dy/dr, 
where r is the dimensionless proper time and c = 1. The particle motion is, then, given by 
the following two equations: 

= ^7 + 2 7*^, (B2a) 
dv y 

= -jp + l\t UV + Yyy U y U y . (B2b) 

The first equation can be solved for u 1 : 

u,= ^--kY- (B3) 



Inserting this u into the second equation gives 



d 2 y , V^^H+^W) UyV _ I 

dT " 2y (y 2 + yl) (-£ + ^yj+tf) \ dr ) 2y (i - ^/yJ+7) 



Finally, replace in flB4j) 

dy 1 dy 2 d 2 y 1 d 2 y 2 1 / dy 2 



dr 2y dr ' dr 2 2y dr 2 4y 3 \ dr 
to obtain 



d 2 y 2 1 L M+v 1 




dr 2 A(y 2 + y 2 )y + ^gq— ~2 ) \dr ) I _ ^yJTf 

which remains finite along the trajectory as long as I < yo- 



(B5) 



(B6) 
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